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ABSTRACT 

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 
0.3/im and 2.4/im. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 
4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative 
photometric calibration is better than 1% across the field of view. The absolute photometric accuracy 
from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy 
spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well 
with the expected values. Finally, as the first step in the scientific analysis of these data we construct 
panchromatic number counts which confirm that both the geometry of the universe and the galaxy 
population are evolving. 

Subject headings: cosmology: observations — galaxies: evolution — cosmology: surveys — cosmology: 
large scale structure of universe 
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1. INTRODUCTION 

Advances in astronomy are often driven by improved 
accuracy and precision along with increases in sensi- 
tivity and area of the available data. The Ca nada- 
France Redshift Survey (CFRSK Lil ly et al.l[l995Th the 
Hawaii Deep Surveys (HDSs) (ICowie et al.1 11999 ) , and 
the Hubble Deep Fie lds (HDFs) (jWilliams et al.l Il99l 
ICasertano et al.l I2000D were the first deep imaging and 
spectroscopic surveys aimed at understanding galaxy for- 
mation and evolution. These discovered the global de- 
cline in star formation at z < 1 and showed that this 
was due to star f ormation occurring in smaller galax - 
ies at later times (|Lillv et al.l Il996l: ICowie et all 119991 ). 
a phenomenon often r eferred to as " Cosmic Downsiz- 
ing" . At the same time iSteidel et all (|1996l 119991 . 12003D 
used the Lyman-Break Galaxy (LBG) color selection 
technique to identify galaxies at high redshift, dramati- 
cally improving the efficiency of spectroscop ic surveys at 
z > 3. Other selections such as the B zK (jDaddi et al.l 
l200l . BX/BM (lAdelberger et al]l200l . and DRG (Dis- 
tant Red Galaxy) (jFranx et al.l 120031 ) have allowed for 
efficient sorting of 1 < z < 3 galaxies. 

Photometric redshifts are the logical extension of color 
selection by estimating redshifts and spectral energy dis- 
tributions (SEDs) from many photometric bands. Un- 
like color selection, photometric redshifts take advantage 
of all available information, enabling redshift estimates 
along with the age, star formation rate (SFR) and mass. 
Unfortunately, photometric redshifts are also suscepti- 
ble to systematics in all bands. This increases the cali- 
bration requirements, especially the required photomet- 
ric accuracy, for modern cosmological surveys such as 
the Great Observatorie s Origins Deep Survey (GOODS) 
(Giavalisco et al. 2004]), the Galaxy Evolution from Mor- 
pholo gy and Spectra l Energy Distributions (GEMS) sur- 
vey (|Rix et al ll2004fh and the Cos mic Evolution Survey 
or COSMOS (|Scoville et al.lf2007bh . 

GOODS and GEMS are designed to study evolution of 
galaxies with look back time, whereas COSMOS is de- 
signed to probe the evolution of galaxies in the context of 
their large-scale structure out to moderate redshift. The 
desire to study large-scale structure in COSMOS neces- 
sitates a 2 square degree area with deep pan-chromatic 
data. Such data have been collected at nearly every ob- 
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servable wavelength from the X-rays to the radio. The 
study of large-scale structures places strong calibration 
requirements on the COSMOS data; for example spatial 
variations in photometry and astrometry must be kept 
to a minimum, typically less than 1% for photometry to 
ensure high quality photometric redshifts and 0.01" for 
astrometry to enable measurements of weak lensing and 
correlation functions. Meeting these calibration require- 
ments is often difficult as multiple instrument pointings 
are used to cover the field. 

This paper concentrates on the ground-based data re- 
duction, the multi-band optical and near-infrared cata- 
log and the steps taken to ensure a high level of pho- 
tometric consistency. The observing strategy for the 
Subaru Suprime-Cam observations, which form the bulk 
of our ground based data, are discussed separately in 
iTaniguchi et alj (|2007l ). In addition, the absolute pho- 
tometric and astro metric system used here is defined in 
lAussel et al.1 (|2007f ). 

An ov erview of the COSMOS project and its goals are 
given in IScoville et al.l (|2007bl ). Details of the Hubble 
Space Telescope (HST) observations, including the Ad- 
vanced Camera for Surveys (ACS), the Wide Field Plane- 
tary Camera 2 (WFPC2), and the Near Infrared Camera 
and Multi-Ob j ect Spe ctrometer (NICMOS) are found in 
IScoville etai] (|2007ah . The ACS data acquisiti on and 
reduction are detailed in iKoekemoer et "all (|2007l ). and a 
monochromatic catalog b ased only on the HST-AC S ob- 
servations is presented in iLeauthaud et all (|2007h . Ob- 
servations at other w avelengths consist of : X-ray ob- 
servations with XMM (Hasingcr et al. 20 0?]), ultraviole t 
(UV) observations with GALEX (jZamoiski et al.ll2007l) . 
mid-in frared observations with the Spitzer Space Tele- 
scope (jSanders et al.l l2007h . sub- mm obs ervations from 
the C altech Sub-mm Observatory (CSO) (|Aguirre et al.l 
|2007| ) and Institut de Radioastronomie Mil limetrique 
(IRAM) 30m telescope (|Bertoldi et al.l 120071 ). and ra- 
dio observations with the Very Large Array (VLA) 
(|Schinnerer et al.ll200l [2007h . 

We begin by presenting an overview of the various data 
sets and photometric systems, the imaging data products 
and the data reduction in Section [2] PSF matching is 
covered in Section [3l and the generation of a multi-color 
catalog is presented in Section |4] Finally, in Section [5] 
we conduct several quality checks, and suggest several 
corrections to the absolute photometry. 

2. OBSERVATIONS AND DATA REDUCTION 

The present COSMOS data were collected on a vari- 
ety of telescopes and instruments, as well as from the 
Sloan Digital Sky Survey (SPSS) s econd data release 
(DR2) archive (|Abazajian et al.l2 004). This paper covers 
the processing of the d ata obtained with Suprime-Cam 
(|Komivama et al.l l2003h on th e Subaru 8.3m telesc ope, 
Megaprime (|Aune et all 120031 : iBoulade et ail 120031 ) on 
the Canad a France Ha waii 3.6m Telescope (CFHT), 
Flamingos (|Elstonlll99cl on the Kitt Peak National Ob- 
servatory 4m ( KPNO-4m), and th e Infrared Side Port 
Imager (ISPI) (jProbst et al.ll2003l ) on the Cerro Tololo 
International Observatory 4m (CTIO-4m) during the 
2004-2005 observing season. 

The telescopes and instruments used for the COSMOS 
survey are presented in Table [TJ A survey efficiency is 
given for each telescope-instrument pair to allow com- 
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Fig. 1. — Filter transmission profiles are shown for the COSMOS 
optical data set from CFHT, Subaru, and HST as of April, 2005. 
These profiles are normalized to a maximum throughput of one 
and include the transmission of the Atmosphere, the telescope, the 
camera optics, the filter, and the detector. The HST F475W data 
only co vers the cent ral 9' X 9', detai l s are given in lScoville et all 
H2007aO . The S PSS lAbazaiian et aTl H2004I) and Johnson-Cousins 
system used by Landolt ( 19921) are shown for comparison. Notice 
the significant differences between the Johnson-Cousins, SDSS, and 
other systems. Color conversions are clearly needed to transform 
from the C OSMOS system to the standard star systems. These 
are given in lAussel et all J2007D . 



parisons between the various data sets. The survey ef- 
ficiency is defined as the telescope collecting area mul- 
tiplied by the detector imaging area (degrees 2 m 2 ) and 
does not include variations in detector sensitivity, sky 
background, or field geometry. This number is most use- 
ful for comparing observations taken in similar bands. 
The filter transmission profiles, including atmospheric 
transmission, telescope reflectivity, instrument optical 
transmission, filter transmission and detector sensitivity 
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Fig. 2. — The transmission profile of the K s band filter from 
KPNO-FLAMINGOS and CTIO-ISPI is shown. This profile is nor- 
malized to a maximum throughput of one and includes the trans- 
mission of the atmosphere, the telescope, the camera optics, the 
filter, and the detector. 
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Fig. 3. — The background limited depth of COSMOS observa- 
tions in the Ultraviolet, Optical, and Infrared are shown. CFHT, 
KPNO/CTIO, and Subaru depths are 5<r in a 3" aperture. SDSS 
depths are those quoted in Abazaiian et al. (2003). The depth of 
the HST- ACS observations is given for a 3" aperture and a 0.15" 
aperture with a point source. A 3" aperture is optimal for color 
measurements, while the 0.15" aperture is the 5c detection limit 
for point sources. The Spitzer-IRAC depths are those expected at 
5<r in a 3" aper ture for observa t ions t aken in 2006. The GALEX 
depths are from Zamojski et al. (2007). 



are plotted in Figures [T] & [2] in units of relative detector 
quantum efficiency normalized to 1 at the peak. 

The Suprime-Cam, Megaprime, SDSS Photometric 
and SDSS Survey cameras have filters distinct from each 
other and from the Landolt standard star system. Even 
the SDSS photometric telescope and SDSS Survey tele- 
scope filter sets differ from one another by 2-4%. To 
differentiate between these filter systems we use a + su- 
perscript for the Suprime-Cam Sloan filters and a * su- 
perscript for the Megaprime Sloan filters; no superscript 
is used for the SDSS survey filters. The designation U, 
B, V, R, and / are used for the Landolt- Johnson-Cousins 
set while Bj and Vj are used for the Suprime-Cam John- 
so n set. Conver s ions b etween these systems are discussed 
in lAussel et~ai1 (|2007f) . 
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TABLE 1 

Telescopes used for COSMOS optical/IR data in 2005-2006 



Telescope 


Telescope 


Instrument 


Field of 


Instrument Wavelength 


Survey 1 


Filters 




Diameter 




View 


Sensitivity 


Efficiency 


used 


CFHT 


3.6m 


Mega-Prime 


56.4' x 57.6' 


3200-llOOOA 


9.19 


u* , i* 


CTIO 


4m 


ISPI 


10.2' x 10.2' 


0.9-2.5/wn 


0.37 


K a 


HST 


2.5m 


ACS-WFC 


3.4' x 3.4' 


4000-11000A 


0.02 




KPNO 


4m 


FLAMINGOS 


10.8' x 10.8' 


0.9-2.5/wn 


0.41 


K a 


SDSS 


2.5m 


SDSS 


25 x 13.5' x 13.5' 


3200-llOOOA 


2.49 


u, g, r, i, z 


Subaru 


8.3m 


Suprime-Cam 


34' x 27' 


4000-llOOOA 


13.8 


Bj, Vj, g+, r+, 














i+, z+, 7VB816 



Defined as the telescope collecting area multiplied by the imaging area in square degrees. 



The wavelength range, depth, and image quality for 
all data presently reduced and included in the version 
2.0 optical/IR catalog are given in Tableland plotted 
in Figure [3J The depth quoted is for a 5c measurement 
in a 3" aperture of an isolated point source at the me- 
dian seeing given in Table [U This should be viewed 
as an optimistic estimate since most objects are ex- 
tended and many are c onfused with neighboring sources. 
iTaniguchi et al.l (|2007t ) present a discussion of detection 
sensitivities and completeness for various Subaru filters. 
The median photometric depths in the COSMOS i + se- 
lected catalog are discussed in Section |4j Table [2] also 
gives a first order offset to the Vega system; however, 
a color term must be applied to get the t rue Landolt- 
Vega s ystem magnitudes-these are given in lAussel et al.l 
(|200l . 

2.1. Data Products 

We took special care in producing data products that 
simplify analysis and are tractable on contemporary com- 
puters 48 . To do this we defined a common grid of sub- 
images for all data products. The starting point for this 
grid is t he COSMOS astrom etric catalog that covers 4 sq 
degrees (|Aussel et al.ll200"7t ) and is larger than all present 
or planned COSMOS data sets. The area is divided into 
144 sections of 10' x 10', each section is covered by an 
image of size 4096 x 4096 pixels with a pixel scale of 0.15" 
. Therefore, adjacent tiles overlap each other by 14.4" on 
all sides. As a result, the vast majority of objects can be 
analyzed on a single image. The layout of the image tiles 
is shown in Figure [H The pixel scale was chosen to be an 
integer multiple of the 0.05" scale used for the HST- ACS 
images. All images and noise maps are scaled to units 
of njy per pixel, which corresponds to a magnitude zero 
point of 31.4. 

For each Subaru and SDSS band an image with the 
original Point Spread Function (PSF) and a PSF ho- 
mogenized across the field within that band is provided. 
For the Subaru Bj, r + , and i + bands, which have excep- 
tional image quality (0.5" - 0.8" seeing), a "best seeing" 
image is also provided. The CFHT images were taken 
in queue observing mode, ensuring a consistent PSF for 
all observations, so only an original PSF image is pro- 
vided for these data. Finally, due to the large variation 
of the PSF in the CTIO and KPNO data, only a PSF 
homogenized image is provided. 

48 All data products discussed in 

this paper are publicly available at 

[j lfittp: / /irs a.ipac. caltech.edu/data/COSMOS/l http: //irsa .ipac.caltech 
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Fig. 4.— The layout of image tiles for the COSMOS field. 



RMS noise maps are also provided for each filter. 
These are on the same tiling scheme and flux scale as the 
images. The RMS maps include noise contributions from 
photon noise, background subtraction, flat fielding, de- 
fect masking, saturation, and cosmic ray removal. They 
do not include the photon noise contribution from object 
flux. 

2.2. Subaru Suprime-Cam 

The Suprime-Cam instrument (|Komivama et al] [2003) 
on the Subaru 8.3m telescope has a 34' x 27' field of view. 
The camera has 10 2k x 4k Lincoln Labs CCD detectors 
which have good sensitivity between 4000A and 10000A. 
Nine Suprime-Cam pointings were required to cover the 
COSMOS field. During 2004 and 2005, data were ob- 
tained in the Bj, Vj, g + , r + , i + , and z + broad-band and 
the NB816 narrow-band filters. These Suprime-Cam ob- 
se rvations, which require d special planning, are detailed 
in lTaniguchi et all (2007). Further observations in eleven 
300A intermediate bands, IA427, IA464, IA484, IA505, 
IA527, IA624, IA679, IA709, IA738, IA767, IA827, and 
one narrow band NB711 were obtained in 2006 and 2007. 
These new observations have been reduced using the pre- 
scription described here, but will be presented elsewhere. 

Objects brighter than 19 th magnitude are saturated 
in a typical exposure and under good seeing the satu- 
ration level can drop to 22 nd magnitude in long expo- 
sures. As a result, it is extremely difficult to astromet- 
ediu^iaily/CXbrlwalS/these data against external astrometric 
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TABLE 2 
Data Quality and Depth 



Filter 


Central 


Filter 


Seeing 


Depth 1 ' 2 


Saturation 2 


Offset froi 


Name 


Wavelength (A) 


Width (A) 


Range (" ) 




Magnitude 


Vega Syst 


a 


3591.3 


550 


1.2-2.0 


22.0 


12.0 


0.921 


a 


o /y / .y 


( ZU 


u.y 


9fi A 

ZD. 4 


10. 


U.ooU 


1 3 

t>J 


A A KCi 7 


oy / 


u.4-u.y 


Z / .O 


1 Q 7 

lo. / 


-U.lol 


9 


4 / Zo. 1 


i inn 


19 17 


99 9 

zz. z 


19 
1Z.U 


117 
-U. 1 1 1 


„4- 


4779.6 


1265 


0.7-2.1 


27.0 


18.2 


-0.117 


v> 


5483.8 


946 


0.5-1.6 


26.6 


18.7 


-0.004 


r 


6213.0 


1200 


1.0-1.7 


22.2 


12.0 


0.142 


r+ 


6295.1 


1382 


0.4-1.0 


26.8 


18.7 


0.125 


i 


7522.5 


1300 


0.9-1.7 


21.3 


12.0 


0.355 


i+ 


7640.8 


1497 


0.4-0.9 


26.2 


20.0* 


0.379 


i* 


7683.6 


1380 


0.94 


24.0 


16.0 


0.380 


F81AW 


8037.2 


1862 


0.12 


24.9+ 


18.7 


0.414 


NB816 


8151.0 


117 


0.4-1.7 


25.7 


16.9 


0.458 


z 


8855.0 


1000 


1-1.7 


20.5 


12.0 


0.538 


z+ 


9036.9 


856 


0.5-1.1 


25.2 


18.7 


0.547 


K s 


21537.2 


3120 


1.3 


21.6 


10.0 


1.852 



1 5ct in a 3" aperture for an isolated point source at the native seeing. 2 In AB magnitudes."^ AB magnitude 
— Vega Magnitude 4- Offset. This offset does not include the color conversions to the Johnsons-Cousins sys- 
tem used by Landolt (1992). Compact objects saturate at i+ < 21.8 due to the exceptional seeing. + The 
sensitivity for photometry of an optimally extracted point source is 27.1, for optimal photometry of a l" 
galaxy it is 26.1 



catalogs such as the SPSS (lAbazaiian et alj 12004 ) and 
USNO-B1.0 (|Monet et all I2003D which only reach 21 s * 
magnitude. To mitigate this limitation, a series of short 
exposures were taken in each band. 

2.2.1. Initial Calibration 

The first step in our Suprime-Cam data reduction is to 
measure a bias level from the over-scan region and sub- 
tract it from all the images. Then all bad or saturated 
pixels are masked. Next, a median bias frame is con- 
structed from over-scan corrected frames. Following the 
over-scan correction, this bias frame is subtracted from 
data and flat frames in order to remove bias structures. 
In particular, the bias level increases near the edges of 
the CCDs farthest from the readout register. 

A median dome flat for each band is then constructed 
from 10-20 bias subtracted flat-field images. The me- 
dian dome flats and median biases are inspected for bad 
pixels, charge traps, and other defects that need to be 
masked. The appropriate median dome flat, with all de- 
fects masked, is then used to normalize all data frames. 
Finally portions of the image vignetted by the guide 
probe are masked from all data frames using the posi- 
tion of the guide probe recorded in the image header. 

After the initial calibration catalogs are generated for 
every data frame with the IMCAT 49 'hfindpeaks' and 
'apphot' routines using 5" diameter apertures. The large 
aperture is chosen to minimize photometric variations 
caused by changes in seeing. This catalog is then used to 
generate an object mask for each frame and to calculate 
the astrometric solution. 

The night sky is subtracted in a two-step process to 
account for fringing and scattered light. First, a nor- 
malized median sky frame is constructed for each night. 
The median frame is generated by masking objects in 
all frames, normalizing every frame to the same median 
flux, and finally median combining the normalized im- 



' [j flrttp: / /www.ifa.hawaii.edu~kaiser/imcat /] http : / / www. if a. hawaii . edu / ~kaiser /imcat / 



ages. The median sky is then scaled to the median back- 
ground level in each frame and subtracted. This removes 
both night sky illumination and fringing. 

After subtracting the median sky residual scattered 
light is visible on the images. This residual light affects 
both the overall flat-field and the background of each in- 
dividual frame. A correction to the flat-field is described 
in Section 12.2.21 while the light in each frame is sub- 
tracted by masking objects and measuring the median of 
the residual background in 128 x 128 pixel squares. A 
background image is then generated by tesselating over 
the grid of medians. After subtracting this background, 
no visible sky structure is left on the individual frames. 
However, this step creates negative halos around bright 
stars and very extended galaxies due to imperfect mask- 
ing. Fortunately, the amplitude of these haloes is similar 
in all frames and can be accounted for as a residual back- 
ground in the combined images. 

After sky subtraction, an astrometric solution is calcu- 
lated separately for all exposures and CCDs by matching 
the object catalogs to the COSMOS astrometric catalog 
(|Aussel et al.1 12007| ) using a 4th order two-dimensional 
polynomial. The polynomial fits are improved by remov- 
ing mismatched objects in an iterative fashion until the 
solution converges (typically in 2 iterations). The result- 
ing scatter between the fit positions and the final astrom- 
etry is always less than 0.2" at the lcr level, independent 
of position. 

Using the astrometric solutions, defects around charge 
bleeds from saturated stars are masked using a list of 
bright stars from the SDSS and the USNO-B1.0. Cos- 
mic rays events are removed by detecting sharp edges 
in the images. Finally, every frame is visually inspected 
to remove internal reflections, satellites, asteroids, and 
other false objects. Once all masking is complete a new 
photometric catalog is generated containing only isolated 
objects in unmasked regions. 

2.2.2. Scattered Light Correction 



() 



Mechanical and optical constraints make it impossible 
to baffle wide field cameras against all scattered light. 
The scattered light is equivalent to an unknown dark 
current added to each image, and must be subtracted 
rather than divided out. As a result the usual flat fielding 
technique of observing a uniform light source such as the 
dome or sky is inaccurate at the 3-5% level. 

For Suprime-Cam the scattered light pattern and 
strength change significantly with the lighting conditions 
and telescope position. Variations as large as ±5% are 
observed at the edges of the field between dark, twilight, 
and dome conditions. Figure [5] shows the difference be- 
tween two dome flats taken at different rotation angles. 
This effect is similar in amplitude and pattern to that 
observed with the 12k and Megacam cameras on the 
CFHT 50 . 

Following the example of CFHT we calculate the true 
flat by observing objects at multiple positions on the 
camera. The true flat can then be solved for as the flat- 
field, which yields the same background subtracted flux 
for an object at any position in the field of view. In 
practice the flat image is generated by dividing the focal 
plane into 128 x 128 pixel regions r, and calculating a 
factor C r for each region. The regions are defined so 
that no region crosses a detector boundary. As a result, 
sensitivity variations due to detectors are also measured. 
We can also allow for an additional factor P e between 
exposures to correct for photometric variations due to 
atmospheric conditions and seeing. However, if the data 
are photometric and corrected for airmass, P e = for all 
exposures. 

For a single object the real magnitude M rea i is de- 
scribed by 



M r 



M + C r + P e 



(1) 



where M is the measured magnitude. If we consider a 
pair of exposures, a and b, we can construct a % 2 relation 
as: 




Fig. 5. — The difference between Suprime-Cam dome flats taken 
in r+ at position angles of and +90 degrees is shown with chip 
to chip sensitivity variations removed. The scale is linear with 
a stretch of -3% to +3% from black to white. Variations in the 
illumination pattern due to scattered light are clearly visible. 




Fig. 6. — The relative correction to the r+ Suprime-Cam dome 
flat is shown with chip to chip sensitivity variations removed. The 
scale is linear with a stretch of -3% to +3% from black to white. A 
correction for scattered light in the vignetted portion of the field 
is clearly visible around the edge of the field of view. 



Nexp Nexp Nobj , 

* 2 = E E E — 

a=0 b=a+l i=o 



Mi.b ~\~ C r a — C r b + Pa 



'Lb 



(2) 

which can be minimized to obtain the C r and P e factors. 
Since an object can only belong to one region in each 
exposure we use the notation C r . a to indicate the region 
an object belongs to in exposure a. 

To ensure that Equation[2]is well constrained a series of 
24 short exposures were taken in each band in each of two 
orientations. Each exposure overlapped its neighbors by 
50%. As a result, any given point in the COSMOS field 
was observed by at least eight different areas of Suprime- 
Cam. These data were calibrated and the photometry 
performed as described in section 12.2. 11 The photometry 
was then corrected for variations in pixel area using the 
Jacobian of the optical distortion. This corrects for the 
assumption that all pixels have the same size when the 
dome flat was applied. Each object is then assigned a 
unique ID by merging the catalogs from each exposure 
with the master astrometry catalog. The IMCAT routine 



'fitmagshifts' was then used to solve Equation [2] for the 
P e and C r factors. 

Two Flat-fields were generated for each band, all the 
data was used in one, and only the data with high photo- 
metric quality in another. No significant difference was 
observed between these flats in any band. As a result 
we included all the data in our solutions. This provides 
an additional check on the photometric calibration be- 
cause an airmass term can be calculated independently 
from the data and from the standard stars. In addition 
non-photometric data can be properly corrected for ex- 
tinction. 



- 1 ['ittp:/ /www. cfht.hawaii.edu/Instruments/Elixir/scattered. html http:/ /www. cfht.hawaii.edu/Instruments/Elixir/scattered. html 



7 



The 24 short exposures were taken while the camera 
was under similar but opposite gravity loads for each ori- 
entation. The flat-fields measured from these two orien- 
tations agree to within 1%, ruling out mechanical flexure 
as a significant source of calibration error. Exposures 
were also taken in the same band over several runs to 
test for changes in the flat-field due to instrument move- 
ments. No significant changes are seen in the flat-field 
between runs spaced as much as a year apart. In con- 
trast, the scattered light pattern, and hence the dome 
flat, changed by as much as 2%. 

Since no variations were observed in the scattered light 
corrected flat-field as a function of telescope position, 
time, or photometric quality, the same flat can be and 
was used for multiple telescope runs. These flats are 
publicly available as part of the COSMOS archive. An 
example of the correction for the r + band is shown in 
Figure [H 

We find that the scattered light component is largest 
in the outer 8' of the field of view which is vignetted 
by the primary mirror. In this outer region the typical 
correction to a dome flattened image is 2-4%. In contrast, 
the central 26' of the field-of-view is stable and flat to 1% 
with no calibration. Therefore the scattered light can be 
safely ignored for the inner regions of Suprime-Cam. 

2.2.3. Photometric Calibration and Image Combination 

After flat fielding object frames taken on photometric 
nights are corrected f or atmospheric ext inction measured 
from standard stars (lAussel et al.ll2007f ). Data taken on 
non-photometric nights are scaled to those taken in pho- 
tometric conditions using the P e factors calculated dur- 
ing the scattered light correction. Exposures with extinc- 
tion greater than 0.5 magnitudes are discarded. Absolute 
photometric calibration is done on the AB system using 
the Subaru filter transmission curves. As a result, all im- 
ages are in units of nJy per p ixel. Color conversi on and 
methodology are discussed in lAussel et al.1 (|2007h . 

After calibration the images in any given band are 
smoothed to the same full-width at half maximum 
(FWHM) using a Gaussian kernel. They are then re- 
sampled onto the final astrometric grid with a sum- 
over-triangles interpolation using the IMCAT 'warpim- 
age' task. Inverse variance-maps, derived from the image 
noise and flat-fields, are also generated and resampled on 
to the final astrometric grid. The variance is scaled so 
that the noise measured in a given sky area on the re- 
sampled images is identical to that measured in the same 
area on the original images. 

Once re-sampled, the images and variance maps are 
combined with the IMCAT 'combineimages' command 
using a weighted sum with outlier pixels, more than 5cr 
from the median, removed. A final RMS map is also 
generated by 'combineimages' that reflects the true pixel- 
pixel RMS. The PSF homogenized images provide the 
most consistent photometry but lose some sensitivity due 
to the smoothing. 

A second combination was done with the original PSF 
images for all bands to provide a maximum sensitivity 
image for detection. Since the PSF varies as a function 
of position and magnitude in these images the aperture 
photometry and colors are less reliable. Finally, for the 
Bj, r + , and i + data a third combination of the images 
was done with only the best seeing data. The resulting 



TABLE 3 
Summary observing log for CFHT 



Telescope Filter Exposure Observation 
Time (h) Date 



CFHT 


u* 


1. 


8 


2003- 


-12- 


-21 






1. 


.8 


2003- 


-12- 


■22 






1. 


9 


2004- 


-01- 


■19 






3. 


4 


2004- 


-01- 


■20 






1. 


8 


2004- 


-04- 


■22 






1. 


.8 


2004- 


-04- 


■25 






0. 


6 


2004- 


-05- 


■22 






0. 


3 


2005- 


-04- 


■04 






0. 


2 


2005- 


-04- 


■05 






1. 


5 


2005- 


-04- 


■09 






1. 





2005- 


-04- 


■11 






1. 





2005- 


-04- 


■14 






1. 





2005- 


-05- 


■04 






0. 


9 


2005- 


-05- 


■05 






0. 


9 


2005- 


-05- 


■06 






0. 


9 


2005- 


-05- 


■08 






0. 


4 


2005- 


-05- 


■29 






0. 


4 


2005- 


-06- 


■02 






0. 


2 


2005- 


-06- 


03 






0. 


4 


2005- 


-06- 


■04 






0. 


4 


2005- 


-06- 


■05 






0. 


2 


2005- 


-06- 


■06 


CFHT 


i* 


0. 


1 


2003- 


-12- 


■21 






0. 


9 


2004- 


-01- 


■15 






1. 


3 


2004- 


-01- 


■17 






2. 


9 


2004- 


■01- 


■18 



images have a FWHM of 0.6" , 0.8" , and 0.5" respec- 
tively. 

2.3. CFHT Megaprime 

The Megaprime camera has a 1 sq degree field of view 
on the 3.6m CFHT telescope. The focal plane is covered 
with 36 2k x 4.5k EEV CCD detect ors with excellent 
response between 3 200A and 9000A (jAune et all 120031 : 
iBoulade et aT]|2003f ). 

Megaprime was used to obtain deep u* band (3798A) 
and shallow i* band images of the COSMOS field. Ob- 
jects as bright as 15 th magnitude are unsaturated in long 
exposures with Megaprime, allowing for an excellent as- 
trometric solution. Observations were taken in sets of 
five dithered exposures forming a five-point dice-face pat- 
tern. A total of five overlapping pointing centers were 
observed to cover the COSMOS field, resulting in data 
four times deeper in the center of the field than on the 
edges. The depth variation is recorded in the noise maps 
discussed in Section 12.11 A summary of the observing 
logs for the Megaprime observations is given in Table [3] 

The CFHT operates in a queue observing mode for 
Megaprime observations which ensures uniform image 
quality and photometry. CFHT and TERAPIX also pro- 
vide a standard reduction pipeline which meets our cali- 
bration requirements. 

Appropriate calibration frames were taken each night 
by the queue observer , and calibrated data were pr ovided 
by the Elixir pipeline (jMagnier fc Cuillandr320 04V This 
pipeline corrects for bias, dark current, flat fielding, and 
scattered light, with the final photometric calibration 
better than 1% across the field of view. 

Further reduction including astrometric and photomet- 
ric calibration, sky subtraction, and image combination 



was provided by the TERAPIX data processing center 51 . 
At TERAPIX the calibrated images provided by Elixir 
were visually inspected with the 'qualityFITS' data qual- 
ity assessment tool and any defective images were re- 
jected. Images with a seeing larger than 1.3" in the i* 
band and 1.4" in the u* band were also rejected. 

All of the images were astrometricall y registered to 
the COSMOS catalog (jAussel et alj|2007l ) using the 'As- 
trometrix' package. They were then resampled using a 
Lancsos-3 interpolation kernel and median-combined us- 
ing the 'SWarp' image combination software. Although 
the median combination is sub-optimal in signal-to-noise, 
it provides the best rejection of cosmic rays. RMS-maps, 
derived from the image noise and flat-fields, were com- 
bined using the same astrometric solution. The image 
scale for the final stack was set to 0.15" per pixel. 

The u* band images were processed twelve months af- 
ter the i* images. The reduction was similar, except a 
new TERAPIX tool 'Scamp' was used to compute the 
global astrometric and photometric solutions for the u* 
images. Finally, a range of quality assessment tests were 
co nducted on the final ima ges similar to those described 
in lMcCracken et all (|2003[ K 

2.4. CTIO ISPI and KPNO FLAMINGOS 

The Flamingos camera on the KPNO 4m |Elstonlll998f) 
and the ISPI camera on the CTIO 4m ([Probst et all 
2003) both provide a field of view slightly larger than 
10' x 10'. These cameras contain a single 2k x 2k 
HAWAII- 2 infrared array detector which is sensitive be- 
tween 0.9/im and 2.4/im. Data from these instruments 
were combined to obtain a K s band image covering the 
entire COSMOS field. However, due to weather and the 
instrument field of view, the depth varies with position. 
The average KPNO exposure time is 1596 sec in K s band 
over the field. With CTIO an exposure time of 1436 
sec in K s was obtained over the whole field. The vari- 
ation in depth is recorded in the noise map discussed 
in Section 12. H A summary of the observing log for the 
KPNO/CTIO observations is given in Table H 

Eighty-one KPNO or CTIO pointings were required 
to cover the entire COSMOS field. Every position was 
covered at least four times with KPNO and three times 
with CTIO. A second grid of 64 pointings offset by half 
a pointing was taken with CTIO to ensure photometric 
consistency. The central nine pointings of the COSMOS 
field were covered with additional passes at the end of 
the 2004 KPNO run. At each pointing a rotated five- 
point dice-face dither pattern with a 1' diameter was 
used. The central position was offset randomly by a few 
arc-seconds between passes to reduce the likelihood of 
bad pixels falling on the same portion of the sky. 

The data were reduced using IRAF 52 with a full 
double-pass reduction algorithm. In the first pass sky 
flats are produced from averaged, sigma clipped, subsets 
of 20-30 dark-subtracted images with similar sky levels. 
A cross-check of different sky flats throughout the night 
shows only low level, large scale, slowly changing gradi- 
ents. A global bad pixel mask is generated from a flat 
used to identify the dead pixels and a pair of dark ex- 
posures used to identify hot pixels. The science data 



TABLE 4 

Summary observing log for CTIO and 
KPNO 

Telescope Filter Exposure Observation 
KPNO 



CTIO 





Time (h) 


Date 




3.6 


2004-02-05 


K s 


4.3 


2004-02-06 


K s 


5.9 


2004-02-07 


K s 


5.0 


2004-02-08 


K s 


5.2 


2004-02-09 


K s 


3.8 


2005-03-31 


K s 


3.5 


2005-04-01 


K s 


5.(5 


2005-04-02 


K s 


0.3 


2004-04-05 


K s 


3.2 


2004-04-06 


K s 


3.8 


2004-04-07 


K s 


3.8 


2004-04-08 


K s 


3.7 


2004-04-09 


K s 


3.2 


2004-04-10 


K s 


2.4 


2005-03-29 


K s 


3.4 


2005-03-30 


K s 


2.8 


2005-03-31 


K s 


2.8 


2005-04-01 


K s 


3.0 


2005-04-02 



[11 
[11 



ttp:/ /terapix. iap.fr/soft/ http:/ /terapix. iap.fr/soft/ 
ttp:/ /iraf.noao.e du/| http:/ /iraf.noao.edu/ 



were then dark subtracted and flat-fielded. Accurate 
positional offsets were determined using IRAF's 'Ima- 
lign' task using multiple well-detected s ources identified 
by SExtractor (|Bertin fc Arnoutsl[l99ll that were com- 
mon to each dithered data set. The images were then 
stacked using integer pixel offsets with IRAF's imcom- 
bine task. These initial stacks of the science data were 
used to generate relatively deep object masks through 
SExtractor's 'CHECKIMAGE.TYPE = OBJECTS' out- 
put file. These object masks were used to explicitly mask 
objects when re-generating the sky flats and in the sec- 
ond pass reduction. Supplementary masks for individual 
images were made to mask out satellites or other bad 
regions not included in the global bad pixel mask on a 
frame-by-frame basis. 

In the full second pass reduction algorithm we individ- 
ually subtracted the sky from each science frame with 8- 
10 temporally adjacent images. Each sky-subtracted im- 
age was then flat-fielded with the corresponding object- 
masked sky flat and any residual variations in the sky 
removed by subtracting a constant to yield a zero mean 
sky level. These individual sky-subtracted images were 
then masked using a combination of the object mask and 
any supplementary mask to remove the real sources and 
bad regions in the sky frames and averaged with sigma- 
clipping to remove cosmic rays. These images were fur- 
ther cleaned of any non-constant residual gradients as 
needed by fitting to the fully masked (object + supple- 
mentary + global bad pixel masks) background on a line- 
by-line basis. The dithered data sets were re-stacked with 
the same offsets determined in the first pass using the 
global and supplementary masks. Finally, an initial as- 
trometric solution was determined using th e 15-50 stars 
on ea ch frame from the USNO-A2 catalog (|Monet et al.l 
EMI). 

After flat-fielding and sky subtraction, frames with 
seeing worse than 1.5" were removed and all remain- 
ing frames were visually inspected. Any frame with es- 
pecially high noise, poor tracking, or other defects was 
removed. Next, every frame was registered to the COS- 
MOS astrometric catalog with a 4th order, 2-dimensional 
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Fig. 7. — Photometric offsets for stars are shown between the SDSS r and Subaru r+ data as a function of RA (left) and DEC (right). 
Notice there are no systematics in RA, but there is a step like pattern in DEC which corresponds to two different SDSS runs. The step 
pattern is due to imperfect PSF matching. 



polynomial. The polynomial fits were improved by re- 
moving mismatched objects in an iterative fashion until 
the solution converged (typically in 4 iterations). The 
resulting scatter between the fit positions and the final 
astrometry was always less than 0.3" independent of po- 
sition. The large astrometric scatter was due to the 0.3" 
pixel size of the detector and poor seeing. 

Each image was then scaled so that its photom- 
etry agrees with th e 2MASS point source catalog 
(jSkrutskie et al.l [2006). A shift of 1.852 magnitudes was 
applied to the 2MASS magnitudes to convert them to the 
AB system. The 2MASS, CTIO, and KPNO K s filters 
are sufficiently similar so that no color terms are needed 
to convert between these filters. After photometric cal- 
ibration, the scatter between overlapping exposures is 
measured to remove any position dependent photomet- 
ric shifts. 

Finally, all data were smoothed with a Gaussian kernel 
to a 1.5" FWHM and combined with 'SWarp', to pro- 
duce images on the COSMOS grid with 0.15" per pixel. 
RMS noise maps, derived from the image noise and flat 
fields, were also combined using the same astrometric so- 
lution. An original PSF image was not produced for the 
K s band data because the PSF variations were too large. 

2.5. Sloan Digital Sky Survey 

SDSS second da ta release (DR2) images 
(|Abazaiian "eTall l200l were used to obtain pho- 
tometry for objects saturated in the Subaru and CFHT 
data. Objects as bright as 10 th magnitude have good 
photometry in the SDSS images. DR2 was used because 
later SDSS data releases contained no new data or 
calibration for the COSMOS field. 

To facilitate photometric measurements, a mosaic of 
the SDSS data was created on the same grid as the other 
COSMOS data. The data were downloaded and cali- 
brated with the "Best" photometric calibration as out- 
lined on the SDSS-DR2 web site 53 . A median sky level 
was then subtracted and a catalog generated for each 
image. The objects in each image were then matched 
to the COSMOS astrometric catalog using a 3 rd order 
2-dimcnsional polynomial with the IMCAT 'fitgeometry' 

53 ['ittp://www. sdss.org/dr2/ http://www.sdss.org/dr2/ 



TABLE 5 

PSF PROPERTIES AND SMOOTHING KERNELS USED FOR 
PHOTOMETRY 



Filter 


PSF FWHM 


Fraction of Flux 


Smoothing Kernel 




(" ) 


in 3" Aperture 


a (" )* 


u 


1.97 


0.716 


0.000 


u* 


0.90 


0.919 


0.662 


Bj 


0.95 


0.942 


0.699 


9 


1.97 


0.725 


0.000 


9+ 


1.58 


0.795 


0.240 


Vj 


1.33 


0.874 


0.521 


r+ 


1.05 


0.914 


0.639 


r 


1.97 


0.708 


0.000 


i 


1.97 


0.709 


0.000 


i+ 


0.95 


0.914 


0.611 


i* 


0.95 


0.891 


0.620 


NB816 


1.51 


0.851 


0.463 


F8UW 


0.07 


0.979 


0.785 


z 


1.97 


0.701 


0.000 


z+ 


1.15 


0.866 


0.585 


Ks 


1.50 


0.759 


0.000 



* FWHM = 1y/2ln(2)cr 



routine. 

The PSF of the SDSS data was homogenized by 
smoothing all images to the same FWHM with a Gaus- 
sian kernel. However, systematic effects of the order of 
5% remain between SDSS 'stripes' due to non-Gaussian 
wings of the SDSS PSF. This is discussed further in Sec- 
tion [3] The COSMOS archive also contains a second 
combination with the un-PSF homogenized images. 

After PSF homogenization, images were re-sampled 
with a linear interpolation using the IMCAT 'warpim- 
age' routine. All images in each band were then com- 
bined with the IMCAT 'combincimages' routine using a 
weighted average. A RMS noise map was also generated 
during the image combination process. The image scale 
for the final stack was set to 0.15" per pixel. 

3. POINT SPREAD FUNCTION (PSF) MATCHING 

A consistent PSF within each band and between bands 
is essential for high-quality photometry. Ideally, all 
bands would have an identical point spread function 
(PSF), but achieving a homogeneous PSF for a data set 
as diverse as COSMOS is extremely difficult due to the 
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TABLE 6 

S Extractor Parameters 



Parameter 


Setting 


Comment 


PARAMETERS-NAME 


cosmos, param 


Fields to be included in output catalog 


FILTER.NAME 


gauss_2 . 5_5x5 .conv 


Filter for detection image 


STARNNW.NAME 


default. nnw 


Neural-Network_Weight table filename 


CATALOG.NAME 


STDOUT 


Output to pipe instead of file 


CATALOG-TYPE 


ASCII 


Output type 


DETECT.TYPE 


CCD 


Detector type 


DETECT.MINAREA 


3 


Minimum number of pixels above threshold 


DETECT .THRESH 


0.6 


Detection Threshold in a 


ANALYSIS-THRESH 


0.6 


Limit for isophotal analysis a 


FILTER 


Y 


Use filtering 


DEBLEND.NTHRESH 


61 


Number of deblending sub-thresholds 


DEBLENDJVHNCONT 


0.0 


Minimum contrast parameter for deblending 


CLEAN 


Y 


Clean spurious detections 


CLEAN.PARAM 


1 


Cleaning efficiency 


MASK.TYPE 


CORRECT 


Correct flux for blended objects 


PHOT_APERTURES 


6.7,13.3,20,26.7,33.3,66.7 


MAG-APER aperture diameter(s) in pixels 


PHOT_AUTOPARAMS 


2.5, 3.5 


MAG-AUTO parameters: <Kron_fact>,<min_radius> 


PHOT.FLUXFRAC 


0.2,0.5,0.8,0.9 


Define n-light radii 


PHOT_AUTOAPERS 


20.0, 20.0 


MAG_AUTO minimum apertures: estimation, photometry 


SATUR.LEVEL 


200000 


Level of saturation 


MAG.ZEROPOINT 


31.4 


Magnitude zero-point 


GAIN 


1 


Gain is 1 for absolute RMS map 


PIXEL.SCALE 


0.1500 


Size of pixel in " 


SEEING.FWHM 


0.95 


Stellar FWHM in " 


BACK.SIZE 


256 


Background mesh in pixels 


BACK_FILTERSIZE 


5 


Background filter 


BACKPHOTO.TYPE 


LOCAL 


Photometry background subtraction type 


BACKPHOTO.THICK 


120 


Thickness of the background LOCAL annulus 


WEIGHT.GAIN 


N 


Gain does not vary with changes in RMS noise 


WEIGHT.TYPE 


MAP.RMS 


Set Weight image type 


MEMORY .PIXSTACK 


5000000 


Number of pixels in stack 


MEMORY J3UFSIZE 


4096 


Number of lines in buffer 


MEMORY.OBJSTACK 


60000 


Size of the buffer containing objects 


VERBOSE.TYPE 


QUIET 


INTERP.MAXXLAG 


3 


Number of bad pixels to interpolate over in X 


INTERP.MAXYLAG 


3 


Number of bad pixels to interpolate over in Y 


INTERP.TYPE 


ALL 


Type of Interpolation 



non-gaussian portion of most PSFs. To ensure a con- 
sistent PSF across the field of view and between bands 
we adopted a two step process. First we homogenize the 
PSF within each band during the data reduction, then 
we match the homogenized PSFs to the band with the 
worst image quality. 

Within each band, we adopt a Gaussian kernel to ho- 
mogenize the PSF between exposures. This works well 
if the seeing variations are small, there are many images 
at the same position to average out the PSF, and the 
photometric apertures are much larger than the seeing 
size. These assumptions, however, break down for large 
seeing variations and small numbers of images. In the 
COSMOS data the effect of non-Gaussian PSF compo- 
nents is negligible (< 0.01 magnitudes) for all but the 
SDSS data when using the PSF matched 3" aperture 
photometry. 

The SDSS data consist of a single exposure at each po- 
sition with seeing variations as large as 1" between the 
two nights which cover the COSMOS field (see Section [2] 
and 12. 5|) . These data are collected in five parallel strips 
in RA with a detector wide gap in declination between 
strips. Two passes are required to completely cover an 
area of the sky. Figure [7] shows the offset between the 
SDSS r and the Subaru r + band photometry measured 
in a 3" aperture as a function of Right Ascension and 
Declination. An offset of 0.06 magnitudes is clearly seen 
in declination between the two SDSS passes due to the 



imperfect PSF matching. This offset is not visible if to- 
tal magnitudes, which correct for seeing variation, are 
compared (see Section I5.1j) . The offsets are similar in 
all SDSS bands, so colors between SDSS bands are not 
significantly impacted. 

For stellar photometry measured with the Subaru data 
in small (~ 1" ) apertures a magnitude-dependent aper- 
ture correction is required between 18th and 21st mag- 
nitude. The correction is due to differences in seeing 
between the short and long exposure data taken with 
Suprime-Cam. The short exposure data were typically 
taken at higher air mass, and hence worse seeing, than 
the long exposure data. In addition, fewer exposures 
were taken with shorter than longer exposure times. As a 
result, the long exposure data are smoothed with a much 
larger kernel than the short exposure data resulting in 
better PSF matching at fainter magnitudes. Corrections 
for th is effect are given for several bands in I Robin et al.l 
(12007ft . which uses a different photometric catalog from 
the one presented here. 

To avoid these PSF matching problems in the multi- 
band catalog, the PSF matching was optimized for a 3" 
aperture. This was achieved by convolving each PSF ho- 
mogenized image with a Gaussian kernel that produced 
the same flux ratio between a 3" and 10" aperture for a 
point source. This method is superior to simply match- 
ing the FWHM since it accounts for the non-Gaussian 
parts of the PSF. However, with this method, only 3" 
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apertures are free from systematic effects. 

We selected the kernel for each image by identifying 
point sources in the ACS images, then using these to 
construct a median PSF for each of the PSF homoge- 
nized images. A Gaussian kernel was then selected which 
yielded the same flux ratio between a 3" and 10" aper- 
ture as the band with the smallest ratio (the K s image) . 
These smoothing kernels are listed in Table [5j 

With the exception of the SDSS images, the worst see- 
ing image was the K s band which contained only 76% of 
the flux in a 3" aperture. To account for the fact that 
the SDSS images have seeing worse than the K s image, 
an aperture correction of -0.06 magnitudes was applied 
to the photometry from these bands. 

4. THE MULTI-BAND CATALOG 

The COSMOS multi-band catalog is derived from a 
combination of the the CFHT i* and Subaru i + original- 
PSF images. The CFHT i* band image alone is too 
shallow, while compact objects in the Subaru i + im- 
age saturate at 21 s * magnitude. Therefore, a combi- 
nation of the two gives the largest possible dynamic 
range for a detection image. The resulting catalog is 
well ma tched in wavelength an d depth to the HST-ACS 
catalog (|Leauthaud et al.ll2007f) and optimally de-blends 
the ground based photometry. The catalog is also op- 
timal for many science goals which require photometric 
redshifts. 

A x 2 image constructed from multiple bands was also 
tried as a detection image but then rejected. The main 
advantage of a x 2 image is pan-chromatic completeness. 
However, the depth and quality of the i + band image 
mean only faint objects with very extreme colors would 
be detected in a \ 2 image but not the i + band image. 
As a result of their faint magnitudes and extreme col- 
ors these objects will have poor photometric redshfits, so 
detecting them is of limited use. Furthermore, the scien- 
tific drivers for a x 2 image strongly favor including the 
K s band which has 1.5" seeing. This poor seeing signif- 
icantly reduces the ability to split close pairs of objects, 
so the benefits of a x 2 image are outweighed by the loss 
in resolution. 

The catalog was generated with SExtractor 
(jBertin fc Arnoutsl I1996T ) run in dual image mode 
on each of the 144 image tiles (see Section 12.11) . Only 
objects in the central 10' x 10' of each tile were recorded 
to avoid duplicate detections in the overlapping regions. 
The combined i + and i* band image was used as the 
detection image for all bands. Absolute RMS maps 
were used as weight maps for both the detection and 
measurement images 54 . Measurements were done on 
the PSF homogenized images, further smoothed to PSF 
match all bands (see Section [3|). 

The number of contiguous pixels, detection thresholds, 
de-blending parameters, and smoothing kernels were ad- 
justed to maximize completeness when visually compared 
to the ACS and x 2 image. Setting these parameters ag- 

54 The default SExtractor settings assume variations in noise 
reflect variations in gain across the image. However, in this case 
both the data and RMS images are in absolute flux units, so the 
gain is 1 everywhere by definition. The SExtractor parameter 
'WEIGHT.GAIN' must be set to N to avoid the default behav- 
ior. Failure to set 'WEIGHT_GAIN=N' will result in incorrect 
error estimates. 



Depth of 



TABLE 7 

Band Catalog Photometry 



r iltcr 


D(T ucpLH 


rtlvlo rvciiigC 


1 

Upper 


2 

Lower 




in 3 /f Aperture 


oi ucpxri 


Quartile 


Quartile 


u 


23.2 


0.3 


22.9 


23.4 


u* 


26.5 


0.2 


26.3 


26.6 


Bj 


26.6 


0.1 


26.6 


26.7 


g 


23.9 


0.2 


23.9 


23.9 


9+ 


26.5 


0.2 


26.5 


26.6 


Vj 


26.5 


0.2 


26.4 


26.5 


r+ 


26.6 


0.2 


26.5 


26.6 


r 


23.6 


0.3 


23.4 


23.8 


i 


22.9 


0.3 


22.8 


23.2 


i+ 


26.1 


0.2 


26.0 


26.2 


i* 


23.5 


0.3 


23.3 


23.7 


NBSW 


25.5 


0.2 


25.4 


25.6 


F8UW 


25.3 


0.1 


25.3 


25.4 


z 


21.5 


0.3 


21.3 


21.7 


z+ 


25.1 


0.2 


25.1 


25.2 


K s 


21.2 


0.3 


21.1 


21.3 



25% of objects with 5<7 measurements are brighter than this 
magnitude^ 25% of objects with 5tr measurements are fainter than 
this magnitude 



gressively results in false detections around bright stars 
and residual image defects. However, these false detec- 
tions can be removed by requiring 5er measurements in 
the detection band, a reasonable FWHM, and a defect 
mask for the detection image. The SExtractor parame- 
ters are given in Table [51 while the software used to gen- 
erate these catalogs is available at the COSMOS website 

55 

Aperture photometry with a 3" diameter is measured 
for each band. This provides the best possible color mea- 
surements by minimizing the effects of PSF variation 
from band to band (see Section[3]). Total magnitudes are 
sub-optimal for color measurements because a correction 
factor must be estimated separately for each band, in- 
creasing the error. Furthermore, SExtractor estimates 
this correction only on the detection image, so estimates 
of total magnitude are only accurate if the image quality 
of the detection and measurement image are identical. 

In contrast, properly PSF matched images have identi- 
cal aperture corrections in all bands, so a single estimate 
of the total magnitude can be used for all bands. We es- 
timate the correction to total magnitude using the offset 
between the total (MAG_AUTO) and the 3" aperture 
magnitude in the detection image. This difference can be 
added to any aperture magnitude to yield a total mag- 
nitude. 

The median 5cr depths for the catalog including PSF 
matching, de-blending, background subtraction, and 
photon noise are given in Table [7] These numbers are 
for total magnitudes and should be used when choosing 
signal to noise cuts and magnitude limits for the COS- 
MOS catalog. An estimate of the RMS variation in the 
5cr limiting magnitude, along with the upper and lower 
quartiles of the 5<r limiting magnitude is also given. 

The general release catalog is cut at a total i + magni- 
tude of 25 and only includes the 2 square degrees with 
uniform multi-band coverage. At this magnitude limit 
and in this area photometric redshifts are reliable, the 

55 [ ittp:/ /cosmos. astro. caltech.edu http://cosmos.astro.caltech.edu 
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catalog is complete, and spurious sources are minimal. 
At fainter magnitudes the catalog begins to be incom- 
plete, have more spurious detections, and photometric 
redshifts begin to behave poorly. A full catalog of all de- 
tections is available on request. However, the full catalog 
should be used with caution. In particular the complete- 
ness and the number of spurious sources will vary as a 
function of position due to differences in the RMS back- 
ground noise. 

4.1. Catalog Contents 

The multi-band catalog contains PSF matched 3" 
aperture photometry and errors for all Subaru and CFHT 
bands along with the KPNO+CTIO K s band and HST 
F814W band data. Non-PSF matched photometry is also 
included for the SDSS bands. Objects with no detection 
are assigned a magnitude of 99 and an error indicating 
the la limiting magnitude expected for that band. Ob- 
jects with no measurement due to lack of coverage, sat- 
uration, or other defects are assigned a magnitude and 
error of -99. 

The catalog photometry uses the photometric zero 
points determined by the standard stars, which are 
known to have systematic offsets on the order of 0.05 
magnitudes and up to 0.2 magnitudes in Bj band. We 
strongly suggest applying the zero point corrections given 
in Table [TBI and discussed in Section IS~3l to get the best 
possible photometry. 

The total (MAG_AUTO) magnitude and the FWHM 
is measured from either the PSF matched Subaru i + or 
the PSF matched CFHT i* band image (if the Subaru 
data is missing or saturated) is included in the catalog 
along with a flag indicating which image was used. As 
discussed in the previous section, since the photometry is 
PSF matched the difference between the total and aper- 
ture magnitudes in the combined i band data provides 
an aperture correction for each object. This aperture 
correction is also included in the catalog. Applying this 
aperture correction to any band will provide a total mag- 
nitude in that band. 

A flag indicating that the object may be a star in- 
stead of a galaxy is also included, however this flag is 
a qualitative assessment and science relying on accurate 
star/galaxy separation should perform a more detailed 
analysis. Quantitative indicators such as the SExtrac- 
tor "CLASS_STAR" parameter are ineffective at sepa- 
rating stars and galaxies because the FWHM varies as 
a function of magnitude on the detection image. Fur- 
thermore, we could not use the HST-ACS data to sep- 
arate point sources because it only covers a fraction of 
the deep ground based data. To overcome this limita- 
tion we manually defined regions on a plot of detection 
image FWHM vs. magnitude which contained the ma- 
jority of ACS point source s (objects with CLASS -STAR 
> 0.9 in the ACS catalog (|Leauthaud et al.ll2007h ). Ob- 
jects falling in this region were then flagged as stars in 
the ground based catalog. 

Flags marking masked regions in the Subaru Bj, Vj, 
i + , and z + images are included in the catalog. These 
masks were generated by visually inspecting images, us- 
ing the SDSS magnitudes as guide to flag bright stars 
and estimate the size of masks. The flag value is the 
area in square arc-seconds of the photometry aperture 
which falls inside a masked region. 



Finally a flag indicating heavily de-blended objects 
is included in the catalog. The default SExtractor de- 
blending settings fail to find objects in areas around 
bright objects due to the high dynamic range of the i + 
band detection image. Unfortunately, de-blcnding ag- 
gressively enough to find faint objects around bright ones 
results in numerous false detections near the bright ob- 
jects. To mitigate this problem objects not detected 
by the default SExtractor de-blending parameters are 
flagged. 

4.2. Catalog Usage Guide 

Objects with all mask and de-blending flags set to 
have the most reliable photometry. However, this re- 
moves a significant fraction of the survey area and some 
flags can be safely ignored under certain circumstances. 

When cross correlating with multi-wavelength detec- 
tions the de-blending and star flags can often be safely 
ignored. However, for clustering analysis based on the 
optical catalog the de-blending flag i s very important, 
especia lly at faint magnitudes (see iMcCracken et all 
(2007b)). Further more, the photometry of faint objects 
de-blended from extended nearby galaxies is suspect due 
to color gradients in the nearby object. 

The photometry masks are specific to the band in 
which they were measured. But since the same instru- 
ment (Suprime-Cam) was used for most of the photom- 
etry they can be safely extended to adjacent bands (for 
instance a combination of the Bj and Vj mask is ap- 
propriate for the g + photometry) . Photometric redshifts 
are affected by masked photometry in a non-linear fash- 
ion, so all photometry masks must be applied to obtain 
a clean photometric redshift sample. However, in spe- 
cific redshift ranges some bands are not as important, so 
the masking could be relaxed if a spectroscopic control 
sample is available. 

4.3. Completeness and Confusion 

For deep surveys the ability to detect objects (com- 
pleteness) and separate superimposed objects (confu- 
sion) are often more important sources of measurement 
error than the photon noise. These quantities are more 
difficult to quantify than the formal noise given in Table 
[2] because they are sensitive to both the data quality and 
the software used. 

The image quality, or seeing has a much larger impact 
on completeness and confusion than the measurement er- 
ror. For a fixed aperture the measurement error increases 
linearly with FWHM, while the peak flux decreases as 
FWHM -1 , and the confusion increases as FWHM 2 . So 
a factor of two difference in seeing corresponds to a 0.3 
magnitude reduction in measurements sensitivity, but a 
0.75 magnitude reduction in peak surface brightness, and 
a factor of four increase in blending. As a result it is of- 
ten possible to achieve much greater flux measurement 
depth than detection depth. SExtractor uses a peak find- 
ing algorithm which is especially sensitive to image qual- 
ity and software settings. The threshold, pixel area, and 
smoothing kernel settings have a large impact on com- 
pleteness while the deblending parameters impact con- 
fusion. These settings are often a compromise between 
detection depth and the number of spurious detections. 

For the COSMOS catalog we chose parameters which 
maximized detection depth but also produced spurious 
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Fig. 8. — The estimated completeness of the i+ band detected 
catalog is shown as a functi o n of magnitude using the method de- 
scribed in McCrackcn ct al. (2003). This method randomly places 
objects with a representative range of morphologies in the field, 
and attempts to recover them. The estimate includes the effects of 
blending and detection completeness. 

detections near the detection limit. These spurious ob- 
jects can be removed with object masks and cuts in mag- 
nitude and FWHM. A similar method to ours is employed 
by the CFHT-LS survey. However, other groups such 
as the Subaru XMM Deep Survey (SXDF), have pre- 
ferred more conservative settings which minimize spuri- 
ous sources at the expense of detection sensitivity. These 
m ore conservative numbe rs are given for all Subaru bands 
in lTaniguchi et all (|2007t ). 

To quantify completeness, simulated objects with a 
representative range of morphologies and magnitudes are 
inserted into the image. SExtractor is then run on the 
image, and the fraction of recovered objects at each mag- 
nitude is measured. No attempt is made to avoid exist- 
ing objects, so the effects of confusion are included in 
the completeness ca lculation. This is ident ical to the 
method described in iMcCracken et al.1 ([2003) . Figure [5] 
shows the results of this simulation for the combined i + 
and i* detection image. The catalog is 91% complete at 
i+ = 25.0, 87% complete at i+ = 26.0 and 50% complete 
at i+ = 27 A. 

5. DISCUSSION 
5.1. Consistency of Photometry 

The COSMOS survey was designed to probe the inter- 
play of large scale structure, and galaxy evolution. High 
quality photometric redshifts are essential for these stud- 
ies. Uncertainties greater than 2 % in redshift will begi n 
to wash out large scale structure (jScoviile et al.ll2007al ). 
so systematic variations in photo-z's need to be less than 
1%. Since a 1% error in photometry typically leads to a 
1% error in photo-z, this places extremely high require- 
ments on the input photometry quality. 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHT-LS) and the SDSS have independently covered 
the COSMOS field, allowing for a quantitative estimate 
of our photometric consistency. The Subaru i + , CFHT 
i*, and HST-ACS F814W data have similar band-passes, 
but were collected on different instruments and reduced 
by different teams using different software, allowing for 
an internal check on photometric consistency. 

The best effort at estimating total magnitude from 




Fig. 9. — The grey scale images show offsets between photometry 
measured by the CFHT-LS and SDSS surveys and COSMOS (left), 
between Subaru i+, CFHT »*, and ACS F814W bands within the 
COSMOS survey (top right), and between 2MASS and COSMOS 
K s band. The images have 1' pixels for all surveys but SDSS 
and 2MASS, which have 2' pixels. The CFHT-LS D2 field only 
covers the central 1 square degree of COSMOS, resulting in a blank 
area. A comparison between SDSS u and CFHT it* is not possible 
because there is no simple linear relationship between these filters. 
The scale is linear in magnitude from black (-0.2 mag) to white 
(+0.2 mag). No variation with position is measurable within the 
measurement errors. 



each survey is used when comparing data to minimize 
effects of aperture size. These specific magnitudes used 
were SExtractor MAG_AUTO values for CFHT-LS, Pet- 
rosian magnitudes for SDSS, and aperture magnitudes 
corrected to total (as described in Section for COS- 
MOS. The CFHT-LS and SDSS photometry were then 
converted to the COSMOS filter system using the con- 
versions given in T ables [5] & [9] and discussed further in 
lAussel et al.l ((2007). Finally the catalogs were merged 
by position with COSMOS allowing for a 1" offset in 
position. 

Figure [5] shows a comparison between the COSMOS 
data and those from the CFHT-LS and SDSS as well 
as between the Subaru i+, CFHT i* , and HST-ACS 
F814W bands for objects brighter than 25th magnitude 
in the COSMOS ground based data and magnitude er- 
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TABLE 8 

Color Conversion Between SDSS and COSMOS Photometric 

Systems 



COSMOS Filter 


Valid for 






Conversion 1 
Equation 




Bj 


-1< (g-r) <1 


Bj 


= 9 


4-0.240(g-r) 


+ 0.029 


9+ 


-1< (g-r) <1 


9+ 


= 9 


- 0.056(g - r) 


- 0.006 


Vj 


-1< (g-r) <1 


Vj 


= 9 


- 0.617(gr -r) 


-0.021 


r+ 


-K (p-r) < 1 


r+ 


= r 


- 0.037(9 -r) 


4- 0.003 


i+ 


-1 < (r - i) < 1 


i+ 


= i - 


- 0.106(r -i) 4 


- 0.007 


z+ 


-1 < (i - «) < 0.8 


z+ 


= z 


-O.llO(i-z) 


f 0.008 



1 A comparison between SDSS u and CFHT u* is not possible because there 



is no simple linear relationship between these filters. 



TABLE 9 

Color Conversion Between CFHT-LS and COSMOS Photometric Systems 



COSMOS Filter 


Valid for 


Conversion 






Equation 



u* 


all 


Same Filter 




Bj 


-0.2 < (9* - r*) < 0.6 


Bj =9* 4-0.432(g* - r*) 


+ 0.047 


g+ 


-1.0 < (9* - r*) < 0.6 
-0.2 < (9* - r*) < 0.6 


9+ = g* 4-0.094(9* - r*) 


4- 0.008 


Vj 


Vj = 9* - 0.545(g* - r*) 


- 0.016 


r+ 


-1.0 < (9* - r*) < 0.9 


r+ = r* - 0.021(9* - r*) 


f 0.001 


i+ 


-1.0 < (9* - r*) < 0.8 


i+ = i* - 0.020(9* - r*) 4 


- 0.005 


F814VK 


0.5 < (i* - z*) < 0.0 


F814W = 2* 4-0.632(i* - 


2*) - 0.116(i* - 2*) 2 - 0.001 


z+ 


-1.0 < (i* - z*) < 0.8 


2+ = z* - 0.128(i* - 2*) - 


- 0.004 



rors smaller than 0.21 (5cr) in the comparison data. A 
comparison between SDSS u and CFHT u* is not possi- 
ble because there is no simple linear relationship between 
these filters. No systematic effects as a function of posi- 
tion are measurable within the photometric uncertainty 
for the 1' bins used for CFHT-LS and 2' bins used for 
SDSS. Smoothing on the size scale of a Suprime-Cam 
pointing yields a typical RMS variation of 0.01 magni- 
tudes across the field between COSMOS, the CFHT-LS, 
and SDSS (see Table [TO]) . While a comparison between 
the COSMOS CFHT i* and Subaru i+ gives 0.007 mag- 
nitudes of scatter (la), and the ACS-F814W data and 
Subaru i + give a dispersion of 0.003 magnitudes. 

Figure [9] also compares the CTIO and KPNO K s band 
photometry to those from 2MASS. Objects with total 
magnitudes of i + > 16 and at least a 5cr detection in 
K s were used in the comparison. No systematic trend is 
measurable with position and the rms variation between 
COSMOS and 2MASS is 0.02 magnitudes on the scale of 
a CTIO/KPNO pointing. At magnitudes brighter than 
i + < 16 the detection image is saturated, resulting in an 
incorrect aperture correction even though the K s band 
photometry is unsaturated. This is manifested as an 
apparent systematic trend with magnitude between the 
COSMOS catalog and 2MASS photometry. However, no 
trend is visible for objects fainter than i + > 16 and with 
measurements made directly on the K s band image. 

The ACS photometry has the best relative calibration 
of all the data sets due to the lack of atmospheric ab- 
sorption. So the excellent agreement between the ACS 
and Subaru i + photometry as a function of position in- 
dicates that the variations seen between COSMOS and 
the CFHT-LS and SDSS are largely due to flat fielding or 
sky subtraction errors in those other surveys. Neverthe- 




F3 1 A'. Al ■■ 

Fig. 10.— The offset between the ACS F814W band photometry 
and the Subaru photometry is s how n as a function of magnitude 
after applying the offsets in Table ff3l The Subaru i+ and 2 + bands 
were used to estimate F814W photometry. The scatter and zero 
point error is within the expected error of the color conversion. 

less, our photometry appears to be constant across the 
COSMOS field to better than 1%, meeting the science 
requirement for large scale structure studies. 

5.2. Galactic Extinction Correction 

The median galactic extinction in the COSMOS field 
is e(B-V)= 0.0195 ± 0.006, which corresponds to an ex- 
tinction of 0.10 ± 0.03 magnitudes in the u* ban d . The 
estimated galactic extinction from ISchlegel et all (1998) 
is provided for each object in the COSMOS catalog. A 
photometric correction for each band can be determined 
from the galactic extinction multiplied by a filter depen- 
dent factor. These factors are given in Table [Tl] for filters 



15 



TABLE 10 
RMS Variation in the Photometric 
Zero Point as a Function of 
Position Between Surveys 



COSMOS Filter CFHT-LS SDSS 



u* 


0.014 


i 


Bj 


0.009 


0.006 


9+ 


0.010 


0.013 


Vj 


0.010 


0.020 


r+ 


0.013 


0.009 


i+ 


0.007 


0.012 


z+ 


0.010 


0.017 



A comparison between SDSS u and 
CFHT u* is not possible because there is 
no simple linear relationship between these 
filters. 



TABLE 11 
Galactic Extinction 
Corrections 



Filter 




Median A A 


u 


4.996724 


0.097436 


u* 


4.690237 


0.091460 


Bj 


4.038605 


0.078753 


9 


3.791856 


0.073941 


9+ 


3.738239 


0.072896 


Vj 


3.147140 


0.061369 


r 


2.649158 


0.051659 


r+ 


2.586050 


0.050428 


i 


1.989881 


0.038803 


i+ 


1.922693 


0.037493 


i* 


1.922912 


0.037497 


F814W 


1.803909 


0.035176 


NB816 


1.744951 


0.034027 


z 


1.467711 


0.028620 


z+ 


1.435914 


0.028000 


K s 


0.340677 


0.006643 



* A A = k A *c(B-V) 



TABLE 12 
Photometric Offsets Measured from 
Other Surveys 



Filter 


Offset to CFHT-LS 


Offset to SDSS 


u* 


0.035 




Bj 


0.125 


0.11 


9+ 


-0.096 


-0.12 


Vj 


0.040 


0.03 


r+ 


-0.080 


-0.07 


i+ 


-0.093 


-0.10 


z+ 


-0.032 


-0.03 



There is no linear relationship between SDSS u 
and CFHT u* 



used on the COSMOS field. These band-pass dependent 
factors are calculated by integrating the filter response 
fu nction against the galac tic extinction curve pr ovided 
bv iBolzonella et all (|2000f ) , originally taken from lAllenl 
41976ft . 

5.3. Absolute Photometric Zero-Point Corrections 

With typical overheads of 15 minutes per standard 
it is extremely difficult to obtain a sufficient number 



TABLE 13 
Photometric 

Offsets 
Calculated 
with 
Spectroscopic 
Redshfits 



Filter 


Offset 


u 


0.0 


u* 


-0.084 


Bj 


0.189 


9 


0.01 


9+ 


-0.090 


Vj 


0.04 


r 


-0.033 


r+ 


-0.040 


i 


-0.037 


i+ 


-0.020 


i* 


-0.005 


NB816 


-0.072 


F814W 


0.000 


z 


-0.037 


z+ 


0.054 


K 3 


-0.097* 



Measured from 
2MASS, not from 
spectroscopic 
redshifts 



of standard stars o n Suprime-Cam. Wit h three to five 
standards per band (|Taniguchi et al.ll2007ft . our standard 
star calibrat i ons ar e accurate to ±0.05 magnitudes (see 
lAussel et al.l ((2007)) . These offsets are larger than de- 
sired for accurate photometric redshifts. 

Comparisons with the CFHT-LS and SDSS yield the 
zero point offsets given in Table [T2j These were esti- 
mated by comparing total magnitudes for point sources 
between 21 st and 24 t/l magnitude for CFHT-LS and 18th 
and 21st magnitude for SDSS converted to the COSMOS- 
AB photometric system. The color conversions between 
th e surveys are given i n Table [5] & [5] and discussed further 
in lAussel et "all (|2007ft . 

The zero point corrections from the two surveys are 
consistent with one another and have an RMS amplitude 
of ±0.06 magnitudes, which is slightly larger than the ex- 
pected error. However, the zero points calculated in this 
way disagree with the ACS F814W photometry by -0.118 
magnitudes and fail to produce photo metric redshi f ts free 
from systematic errors. Furthermore. Illbert et all (2006) 
find that the CFHT-LS zero-points are inaccurate at the 
0.05 magnitude level. 

Considering the number of present and ongoing obser- 
vations of COSMOS, obtaining spectrophotometric stan- 
dards in the field is a reasonable lon g term solution t o 
obtaining better quality photometry (|Wolf et alj 2001). 
In the interim we are forced to rely on the existing cali- 
brations and spectra of galaxies to re-calibrate the pho- 
tometric zero points for photometric redshifts. 

These zero point offsets are calculated by fitting the 
measured photometry, corrected for galactic extinction, 
to galaxy templates at the known redshift and calculating 
the offset between the measured and expected photom- 
etry. To avoid systematic variations due to calibration 
errors in the galaxy templates, offsets between bands are 
calculated in rest frame wavelength bins of 100A sepa- 
rately for each template, then combined in a weighted 
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Fig. 11. — Four Color-Colo r plots are sho wn for point sources in the COSMOS field (black) a l ong with the expected colors from the 
BaScl 3.1 (Wcstcra ct al. 2002) and NextGen (Hauschildt & Baron 2005; Brott & Hauschildt 2005; Hauschildt ct al. 2002) stellar libraries 
at the expected median metallicity of [Fe/H1=-0.4 (See lRobin et al.1 11 20071 )) (green), and the lPicklesi 1119981) spectral library ([Fe/H]~ 0) 
(red). Notice that the Pickles library predicts systematically redder colors than the observations in the UV due to the high metallicity of 
these stars. 



average. This is effective because calibration errors in 
the template will create offsets which vary as a function 
of rest wavelength in the same way for all bands, while 
zero point offsets between bands will be constant as a 
function of wavelength. Since these offsets are relative, 
the F814W is used as the reference for the absolute zero 
point. The offsets calculated using this method are given 
in Table [13] 

This method is effective with a large sample of spectro- 
scopic redshifts covering a large range in redshift along 
with photometry taken in many adjacent bands. For 
COSMOS there are insufficient data to correct the K s 
band due to the lack of photometric data between 0.9/im 
and 2.2^m. As a result we rely on the 2MASS zero point 
for those data. Although the individual exposures were 
tied to 2MASS, an offset of -0.097 magnitudes is mea- 
sured between the catalog photometry and the 2MASS 
catalog. This offset is likely an aperture correction be- 
tween the photometry measured on the individual CTIO 
and KPNO exposures and the 2MASS catalog. 

With the exception of the Suprime-Ca m Bj band, 
which has known calibration problems (see lAussel et al.l 
(|2007h ). the offsets are within the expected error of ±0.05 
magnitudes. After applying these offsets no systematic 
trend is measurable between photometric and spectro- 
scopic redshifts for 842 objects between < z < 1.2. In 
addition, the i + and z + band agree with the F814W pho- 
tometry to 0.007 magnitudes and the offsets between the 
COSMOS and CF HT-LS photometry agree with those in 
lllbert et al.l (|2006[ ) after applying color corrections. Fur- 
thermore, after applying the offsets the colors of stars 



agree with the predicted colors extremely well (see Sec- 
tion [53J. These tests indicate the zero points are within 
0.01 magnitude of the true AB zero points after applying 
these offsets. 

These corrections were not applied to the released cat- 
alog since they can not be verified with external calibra- 
tion sources at this time. However, we recommend apply- 
ing the photometric redshift offsets for the best possible 
photometry and colors. 

5.4. Star Colors 

Stars are a good diagnostic of color accuracy because 
they form a tight sequence in most optical and near-IR 
color-color plots. Offsets as small as a few hundredths 
of a magnitude are visible when comparing expected and 
measured star colors. Furthermore, star colors are sensi- 
tive to the filter transmission profiles, prov iding a valu- 
able c heck of the instrumental performance ( jAussel et al.l 
l2007h . 

Even at the resolution of HST, compact galaxies and 
quasars contaminate a star selection based on morphol- 
ogy. These objects create scatter in color-color plots, 
obfuscating the stellar locus. The BzK col or-color dia- 
gram provides a much cleaner star selection (jDaddi et al.l 
2004), and was used to select the objects plotted in Fig- 
ure [TT1 The BzK method is biased against faint blue 
stars due to the shallow K s band data, however, the ef- 
fect is minimal since only objects with greater than 10<7 
detections are plotted. 

Stellar libraries e.g. (Pickles 199§) typically con- 
tain solar metallicity stars ([Fe/H]~0), which differ in 
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Fig. 12. — I band number counts are shown for the COS- 
MOS g round based cata l og, th e COSMOS F814W weak lensing 
catalog Leauthaud ct al. (2007), the Hawaii-Hubblc deep field (H- 
HD F-NnCapak et al.1120041), the Hubble Deep Field North (HDF- 
N) nWilliams et a l. 1996; Metcalfe et al. 2001]), the Hubble Deep 
Field South (HDF-S) HMetcalfe et al J 120011) , the Herschel Deep 
Field (Metcalfe ct al. 2001), the S PSS HYasuda et al I l200lT) , the 
Canada France Deep Fields (CFDF1 [McCracken et*aUll25oa) . and 
the Canada-France-Hawa ii Telescope Legacy Survey (CFHT-LS) 
IIMcCra ckcn ct al. 2007b]). Our counts are in good agreement with 
other surveys up to our 80% completeness limit at i + = 26.5. 

color from the dominant sub-solar metallicity ([Fe/H]~- 
0.4), thick disk po pulation in the COSMOS field (see 
iRobin et al.l ([2007)). This results in small offsets, espe- 
cially in the ultraviolet where metal line absorption will 
cause higher metallicity stars to appear redder. 

Figure QT] shows four differen t color- color plots for 
stars along with colors for the iPicklesI (119981) library 
and a combination of the BaSel 3.1 (jWestera et al"1l2002D 
and NextGen theoretical librar i es jH auschild t fc Baron! 
12001 IBrott fc Hauschildtl [20051 : IHauschildt et al.ll2002f) 
at [Fe/H]=-0.4. The star colors agree extremely well in 
all color-color plots with an offset between the median 
expected and actual colors of less than 0.02 magnitudes 
in most bands. A correction of -0.05 magnitudes is in- 
dicated by the stellar track for u* , and systematic dif- 
ferences between the predicted and actual colors for Bj 
and z + are observed. These are likely du e to our lim- 
ited m easurements of the filter throughput (lAussel et al.l 
120071) . 

5.5. Number Counts 

Differential number counts of galaxies are a simple 
but powerful measurement of the geometry of space, the 
evolution of the galaxy population and the evolution of 
structure in the universe. They also provide a valuable 
check on data quality because they are sensitive to pho- 
tometric calibration errors, detection completeness, and 
spurious detections. Figure [T2l shows the I band number 
counts from the COSMOS catalog compared to other sur- 
veys. The numbers are in good agreement, however the 
COSMOS i+ band counts are higher than the F814W 
ACS counts at magnitudes fainter than 24. This is likely 
due to the lack of sensitivity at low surface brightness in 
ACS compared to Subaru and some contamination from 
spurious sources in the ground based catalog. The i + 
band galaxy number counts are given in Table [Til stars 
have been removed using a cut in ground based FWHM 
since ACS does not cover the entire 2 square degrees. 
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Fig. 13. — Number counts from the i + detected catalog are shown 
for u* , Bj, g + , Vj, r+, i + , F814W, z + , and K 3 bands divided by 
a Euclidian slope. The u* counts are nearly flat while K s band 
counts are very steep. A combination of a f2 m = 0.3, Q v = 0.7 
geometry, galaxy evolution, and redshifting is responsible for the 
band dependent behavior. 

Figure [13] show the number counts for the COSMOS 
bands normalized to a slope of 0.6, which is expected in 
a static Euclidian universe. In a universe with no galaxy 
evolution and a cosmology of fl m = 0.3, fl v — 0.7, the 
number count slope should be below 0.6 in all bands. 
Furthermore, galaxies have a flat to red spectrum in units 
of F„, so the k-correction is positive in AB magnitudes. 
This results in galaxies fading faster with redshift than 
they would from simple luminosity distance, further flat- 
tening the number count slope. 

The slope of the normalized number counts in Figure 
1131 steepen with increasing wavelength indicating both 
the universe and the galaxy population is evolving. The 
UV luminosity of a typical galaxy increases with redshift 
due t o increased star formation activity (|Cowie et al.l 
[19991 iLillv et~aTI Il99fl iDahlen et alll2T)0~5h . This galaxy 
evolution counters the effects of an evolving universe, re- 
sulting in a u* band slope clo se to 0.6. In the n ear-IR, 
there is little galaxy evolution (Dahlen et al. 2005) so the 
counts in the K s band are far below a slope of 0.6. 

6. CONCLUSIONS 

Deep ground based data were presented for the COS- 
MOS field in 15 bands between 0.3^m-2.2/zm along with 
an i band selected catalog of objects. We show that 
these data have an extremely high level of photometric 
consistency necessary for scientific pursuits such as large 
scale structure studies. The expected zero-point varia- 
tions are < 1% across the field; this level of photometry 
was achieved by constructing flat fields directly from ob- 
ject fluxes rather than using sky or dome flats. 

Variations in the point spread function (PSF) are a 
significant source of uncertainty in ground based pho- 
tometry. This is due to the fact that the non-Gaussian 
portion of the PSF is difficult to match across multiple 
bands, resulting in 2-5% errors in color measurement if 
uncorrected. The variation was minimized by adjusting 
the Gaussian filter used for PSF matching to provide the 
same fraction of light to fall in a 3" aperture in all bands. 
However, this means that only the 3" aperture photom- 



18 



TABLE 14 
i+ Band Number Counts 
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etry is reliable for colors. Photometry measured in other 
aperture sizes should be corrected for variations in the 
PSF. 

An i + band selection was used for the multi-band cat- 
alog after analyzing the tradeoffs between a \ 2 an d i 
band detected catalog. The decrease in resolution and in- 
creased problems with a x 2 catalog outweigh the benefits 
of increased sensitivity and pan-chromatic completeness. 
The resulting catalog is 91% complete at i + = 25.0, 87% 
complete at i + — 26.0 and 50% complete at i + — 21 A. 



Abazajian et al., 2004, AJ, 128, 502 
Adelberger et al., 2004, ApJ, 607, 226 
Aguirre et al., 2007, submitted to ApJ 

Allen, C. W. 1976, Astrophysical Quantities (University of London, 

The Athlone Press) 
Aune et al., L. G. 2003, Proc. of the SPIE, Volume 4841, pp. 513- 

524 (2003), 513-524 
Aussel et al., 2007, in Prep 
Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 
Bertoldi et al., 2007, submitted to ApJ 

Bolzonella, M., Miralles, J.-M., & Pello, R. 2000, A&A, 363, 476 
Boulade et al., L. G. 2003, Proc. of the SPIE, Volume 4841, pp. 
72-81 (2003), ed. M. Iye & A. F. M. Moorwood, 72-81 



Ou r photometric zero -points measured from standard 
stars (|Aussel et al.ll2007l ) are accurate to 5%. These zero- 
point offsets can be significantly improved by using galax- 
ies with known redshifts to adjust the zero-point calibra- 
tion. The corrected zero-points appear accurate to better 
than 2% based on star colors; this will be verified with 
future internal spectrophotometric standards. 

Number counts measured from the i + band selected 
catalog agree well with previous studies. The effects of 
galaxy evolution and an expanding universe are clearly 
visible in the pan-chromatic counts, demonstrating the 
importance of multi-color surveys. 

We would like to thank the COSMOS team 
http://cosmos.astro.caltech.edu/, the staff at Cal- 
tech, CFHT, CTIO, KPNO, NAOJ, STSCI, Terapix, and 
the University of Hawaii for supporting this work and 
making it possible. Support for this work was provided 
by NASA grant HST-GO-09822 and NSF grant OISE- 
0456439. 

FLAMINGOS was designed and constructed by the IR 
instrumentation group (PI: R. Elston) at the University 
of Florida, Department of Astronomy, with support from 
NSF grant AST97-31180 and Kitt Peak National Obser- 
vatory" 

Funding for the Sloan Digital Sky Survey (SDSS) has 
been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Aeronautics and 
Space Administration, the National Science Foundation, 
the U.S. Department of Energy, the Japanese Monbuk- 
agakusho, and the Max Planck Society. The SDSS Web 
site is |http: / /www.sdss.org/| 

The SDSE is managed by the Astrophysical Re- 
search Consortium (ARC) for the Participating Institu- 
tions. The Participating Institutions are The Univer- 
sity of Chicago, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, The Johns Hop- 
kins University, Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max- 
Planck- Institute for Astrophysics (MPA), New Mexico 
State University, University of Pittsburgh, Princeton 
University, the United States Naval Observatory, and the 
University of Washington. 

D. Thompson acknowledges support from NASA 
through a Long Term Space Astrophysics grant (NAG5- 
10955, NRA-00-01-LTSA-064). 

This work is based (in part) on data products pro- 
duced at TERAPIX data center located at the Institut 
d'Astrophysique de Paris. 



Brott, I. k. Hauschildt, P. H. 2005, in ESA SP-576, 565-+ 
Capak et al., 2004, AJ, 127, 180 
Casertano et al., 2000, AJ, 120, 2747 

Cowie, L. L., Songaila, A., & Barger, A. J. 1999, AJ, 118, 603 
Daddi et al., 2004, ApJ, 617, 746 
Dahlen et al., 2005, ApJ, 631, 126 

Elston, R. 1998, in Proc. of the SPIE Vol. 3354, p. 404-413, 404-413 
Franx et al., 2003, ApJ, 587, L79 
Giavalisco et al., 2004, ApJ, 600, L93 
Hasinger et al., 2007, |astro-ph/06123lT1 

Hauschildt, P. & Baron, E. 2005, Memorie della Societa 
Astronomica Italiana Supplement, 7, 140 



19 



Hauschildt, P. H., Allard, F., Schweitzer, A., & Baron, E. 2002, in 

ASP Conf. Ser. 274, 95-+ 
Ilbcrt et al., 2006, ArXiv Astrophysics e-prints 
Koekemoer et al., 2007, astro-ph/0703095 

Komiyama et al., 2003, Proc. of the SP1E, Volume 4841, pp. 152- 

159 (2003)., 152-159 
Landolt, A. U. 1992, A J, 104, 340 
Leauthaud et al., 2007, |astro^p h/0702359 
Lilly et al., 1995, ApJ, 455, 50 
Lilly et al., 1996, ApJ, 460, L1+ 

Magnier, E. A. & Cuillandre, J.-C. 2004, PASP, 116, 449 

McCracken et al., 2007, astro-ph/0704.2545 

McCracken et al., 2003, A&A, 410, 17 

Metcalfe et al., 2001, MNRAS, 323, 795 

Monet et al., 1998, VizieR Online Data Catalog, 1252, 

Monet et al., 2003, AJ, 125, 984 

Pickles, A. J. 1998, PASP, 110, 863 

Probst et al., 2003, Proc. of the SPIE, Volume 4841, pp. 411-419 

(2003)., 411-419 
Rix et al., 2004, ApJS, 152, 163 



Robin et al., 2007, |astro-ph/0612349] 
Sanders et al., 2007, astro-ph/0701318 
Schinnerer et al., 2004 AJ, 1L>8, 1974 
Schinnerer et al., 2007, |astro-ph/0612314l 

Schlegel, D. J., Finkb einer, D. P., fe Uay ls, M. 1998, ApJ, 500, 525 

Scoville et al., 2007a, astro-ph/0612306 

Scoville et al., 2007b, astro-ph/0612305 

Skrutskie et al., 2006, AJ, 131, 1163 

Steidel et al.,1999, ApJ, 519, 1 

Steidel et al., 2003, ApJ, 592, 728 

Steidel et al., 1996, A J 112, 352 

Taniguchi et al., 2007, astro-ph/0612295 

Westera et al., 2002, A&A, 381, 524 

Williams et al., 1996, AJ, 112, 1335 

Wolf et al., 2001, A&A, 365, 681 

Yasuda et al., 2001, A J, 122, 1104 

Zamojski et al., 2007, |astro-ph/0701478"1 



